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I— i. Abstract 
CLh ■ 

' We study the motion of the steady compressible heat conducting viscous fluid in 

• . a bounded three dimensional domain governed by the compressible Navier-Stokes- 

'5 I Fourier system. Our main result is the existence of a weak solution to these equations 

^ ' for arbitrarily large data. A key element of the proof is a special approximation of 
, Q the original system guaranteeing pointwise uniform boundedness of the density. 

! Therefore the passage to the limit omits tedious technical tricks required by the 

I standard theory. Basic estimates on the solutions are possible to obtain by a suitable 

^ ■ choice of physically reasonable boundary conditions. 

^ : 
^ : 

^ '• 1 Introduction 

; 

I We consider the following system of partial differential equations describing the steady 

O ■ flow of a compressible heat conducting Newtonian fluid in a bounded three dimensional 

>• ■ domain Q 

^ ■ (1.1) diY{gv) = 0, 



(1.2) div{gv (g) -u) - div S{v) + Vp{g, 9) = gF, 

(1.3) div {ge{g, e)v) - div (^(0) V^) = S{v) : Vv - p{g, 9) div v, 

where g : Vt ^ is the density of the fluid, v : i7 — is the velocity field, S{v) = 
2fiD{v) + A(divf)J is the viscous part of the stress tensor, D{v) = |(Vi' + (Vv)'^) is 
the symmetric part of the velocity gradient, p{-, ■) : R([ X M+ ^ ]R+, a given function, 
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is the pressure, F : fl ^ is the external force , e(-, ■) : Mj" x — > M^, a given 
function, is the internal energy. The system (1.1)-(1.3) is known as the compressible 
Navier-Stokes-Fourier equations or the full Navier-Stokes system [6]. 
We assume that the constitutive equation has the form 

(1.4) p{q,6) = aiQ^ + a2Q9, 01,02 >0, 

i.e. the pressure has one part corresponding to the ideal fluid and a so called elastic part; 
for more information see e.g. [6]. Even though we could consider more general pressure 
laws, we restrict ourselves to this simple model to avoid unnecessary technicalities in the 
proof. The corresponding internal energy takes the form 

7— 1 

(1.5) e(^,e) = a2^ + ai^— , 

7-1 

see e.g. [6] or [1]. Note that in the full generality, the equation (1.3) should be replaced by 
the conservation of the total energy, instead of conservation of the internal energy only. 
For sufficiently regular class of solutions, including that we are going to construct, the 
balance of the kinetic energy is just a consequence of the momentum equation. We further 
simplify (1.3). As our solution will be such that q e L^{Q) and v e (Q), p < 00, we 
get due to the fact that div(^i;) = in the weak sense (see [16]) 

div ( — ^ — Q^v) = —g^ diyv, 
^7 — 1 ^ 

again in the weak sense. Thus we may write instead of (1.3) (we put ai = 02 = 1) the 
energy equation in the form 

(1.6) div (gOv) - div (k(^) V^) = S{v) -.Vv-gO div v. 

The viscosity coefficients are for the sake of simplicity considered to be constant such that 
the conditions of the thermodynamical stability 

2 

(1.7) ii>0, A + -/x>0 

are satisfied. Finally, the heat conductivity is assumed to be temperature dependent, i.e. 

(1.8) K{e) ^asil + e""), a3,m>0. 

This fact is important for our study, we are not able to consider constant heat conductivity. 
Our domain Q is sufficiently smooth, at least a domain. Wc supplement the system 
(1.1), (1.2) and (1.6) with the following boundary conditions at dfl. For the velocity, we 
consider the slip boundary conditions 

(1.9) V • n = 0, Tfe • (T(p, v)n) + fvTk^O at dfl, 

where r^, k = 1,2 are two perpendicular tangent vectors to dfl, n is the outer normal 
vector and T{p, v) = —pI+S{v) is the stress tensor. The slip coefficient / is non-negative 
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(if / = we assume additionally that fl is not rotationally symmetric) . Recall that / = 
corresponds to the perfect shp while / — > oo leads to the homogeneous Dirichlet boundary 
conditions. 

Concerning the temperature, we assume that 

f)f) 

(1.10) K(e)— + L{e){e-eo)^o at an, 

where '■ dft — > R"*" is a strictly positive sufficiently smooth given function, say 9o € 
C^idQ), <0^<0o<9* <oo with 9,, 9* e M+ and 

(1.11) L{9) = a4{l + 9^), leR^. 
We must also add the prescribed mass of the gas 

(1.12) / gdx^M>0. 

The objective of this paper is to prove the existence of weak solutions to problem 
(1.1)-(1.12) for arbitrarily large data. Till now only partial results have been proved 
(see e.g. [2], [9], [14], [15]) and only known general theorems concern weak solutions to 
the evolutionary version of the system [6] . The main obstacle was to construct suitable a 
priori estimates. Due to properties of the boundary conditions (1.10) wc are able to obtain 
a nontrivial energy bound for weak solutions, saving the thcrmodynamical structure of 
the system. In the case of the barotropic gas we do not meet such difficulties. The energy 
bound follows elementary from the momentum equation. Unfortunately, it is not the only 
difference. The standard methods introduced by P.L. Lions [9] do not work successfully 
for the heat conducting case. However, a generalization of the technique introduced in 
[11], [17] gives us sufficient tools to solve the stated problem. 

An approach to system (1.1)-(1.12) was considered in the book [9], unfortunately, 
this result can be viewed as conditional only, since instead of (1.12) the author assumed 
artificially that weak solutions satisfy gPdx — for sufficiently large p. On the 
one hand, this condition is physically not acceptable, on the other hand, it simplifies 
considerably the mathematical analysis. Nevertheless, this result shows us what is the 
difference in techniques for the barotropic and heat conducting models. 

Looking on results concerning the classical solutions for problems with small data, we 
realize that the heat conducting system has the same mathematical structure (difficulties) 
as the barotropic version of the model. Thus results from [2], [15] are almost immediately 
transformed to the case of the system (1.1)-(1.12). For large data solutions the energy 
equation starts to play an important role, essentially changing the properties of the whole 
system. 

The evolutionary case of the system (1.1)-(1.12), under general assumptions on the 
pressure law was considered in [8] and [7]; however, the presented technique treats only 
the situation when the fluid is thermically isolated, i.e. = at the boundary. It 
guarantees immediately the energy bound for weak solutions, but considering the limit 
t ^ oo, the only solution which can be obtained as the limit for large times (with time 
independent force) is the solution with constant temperature. This is connected to the 
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fact that the model does not allow the heat transfer through the boundary and either the 
energy increases to infinity (non potential force) or the temperature approaches a constant 
value (potential force). Boundary condition (1.10) allows the heat transfer through the 
boundary, guaranteeing the balance of the total energy, and thus we are able to prove 
existence of solutions which are definitely nontrivial. 
The main result of this paper is the following. 

Theorem 1 Let fl E be a bounded domain in M^. Let F e L^{Q) and 

37-1 

7>3, m^l + l>-^ — -. 

87 — 7 

Then there exists a weak solution to (1.1)-(1.12) such that 

Q e L^{n), V e Wl{Q) and 9 e Wl{n) for alll<q< 00. 

The solution constructed by Theorem 1 is meant in the following sense. 

Definition 1 The triple (g.v.O) is a weak solution to (1.1)-(1.12), if Q E Ls{VL), s > 2'y, 
V e (^); ^ ^ (^) (^nd O'^VO e Li{Vt), v ■ n = at dQ in the sense of traces and 

(1.13) / QvVrj^O Vr;eC°°(n), 

Jn 

/ ( — (g) V : Wcp + 2nD{v) : D{cp) + Adiv vdiv cp — p{g, 9) div (p)dx 

(1.14) J^. . 

+f {v Q t) ■ {ip Q T)da ^ qF ■ ipdx y^p e C'"{fiy,ip ■ n ^ at dfl 
JdQ Jn 

(we denoted by v Q t the vector v — {v ■ n)n) and finally 

[ {k{9)V9 -Vip - q9v ■Vip)dx+ ! L{9){9 - 9o)^pda 

(1.15) p -'^^ 

= / {2n\D{v)\^'il;dx + X{dwvf'ijj- Q9divvjl;)dx V^' G C~(n). 
Jn 

The proof of Theorem 1 will be based on a special approximation procedure described 
in the next section which is the kernel of our method. This section includes also a priori 
estimates for the approximation. The structure of the approximative system gives us 
immediately the approximative density bounded uniformly in Loo, b^it wc must prove 
refined L^o estimates to verify that the limit solves the original system (1.1)-(1.3). This 
idea has already been successfully applied in [11] and [17] in the case of barotropic fiows. 

The third section contains a detailed proof of existence to the approximative system. 
Here the main difficulty comes from the energy equation, since the required positiveness 
of the temperature does not follow immediately. In the next section we introduce an 
important quantity, the effective viscous fiux and prove its main properties, i.e. the com- 
pactness. This feature allows to improve information about the convergence of the density, 
which is the basic/fundamental fact in the theory of the compressible Navier-Stokes equa- 
tions. The last section describes the refined L^o estimates for the approximative density 
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and the passage to the hmit. Then we prove that the hmit is indeed our sought solution 
in the meaning of Definition 1. 

As the reader may easily check, our method works for slightly larger class of the 
pressure laws. It allows to consider e.g. 



where Pb{[)) is a strictly monotone function which behaves for large values as g'^ . The main 
steps of this generalization are similar to the barotropic case and can be found in [17]; 
since our problem is technically enough complicated, we shall avoid such generalizations. 

Our new result is closely related to the barotropic version of the system (1.1)-(1.12). 
Let us remind the state of the art in this theory. The steady compressible Navier-Stokes 
equations for arbitrarily large data were firstly successfully studied in the book [9], where, 
in the case of p{g) = the existence of renormalizcd weak solutions was shown for 7 > 1 
(A^ = 2) and 7 > | {N = 3) for Dirichlet boundary conditions. For potential forces 
with a small non potential perturbation the existence was improved in [13] for 7 > | 
{N — 3). In the recent paper [5] the authors proved the existence in two space dimensions 
also for 7 = 1. Sec also [3], where the authors considered the three dimensional case 
and got existence for certain 7-s less than |, however, for periodic boundary conditions. 
P.L. Lions also considered the existence of solutions with locally bounded density: for 
the case of Dirichlet boundary conditions he was able to show their existence for 7 > 1 
{N — 2) and 7 > 3 (A^ = 3). Nevertheless, to prove Theorem 1 the above methods are 
not sufficient, thus we present our new approach for the heat conducting model. 

Throughout the paper we use the standard notations for the Lcbcsgue, Sobolev, etc. 
spaces; generic constants are denoted by C and sequences e — > always mean suitable 
chosen subsequences ^ O"*". For the sake of simplicity we put ai — a2 — — — 1. 

2 Approximation 

This section contains one of the main difficulties in the proof of Theorem 1 — to find 
a good approximation of the problem (1.1)-(1.12) for which we arc able to show existence 
and prove the corresponding a priori estimates. We present the approximative system as 
well as the proof of the fundamental a priori estimates. Next section deals then with the 
solvability of this system as well as with further a priori bounds. 

Our approximative system will contain two parameters: a number e > and an 
auxiliary function K{-) defined by a number k > as follows: 



moreover we assume that K'{t) < for t e {k,k + 1), where k e W^. In the last section 
we pass with e — > 0"^ and we shall show that we may take k sufficiently large such that 



p{g,9) = pb{g) + g9, 



(2.1) 
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K[g) = 1 for our solution. The approximation of our problem (1.1)-(1.12) reads as follows 

€g + div{K{g)gv) — eAg — ehK{g) 



(2.2) 



where 



^ diy{K{g)gv ® v) + ^K{g)gv ■ Vv - div S{v) + VP{g, 9) = gK{g)F 

-div + +div (v J' K{t)dty + div (^K{g)gvy 

+K(g)gv ■ Ve - eK{g)v ■ Vg = S(v) : Vv ^ 



> in Q,, 



(2.3) P{q,0)= f\t"'-^K{t)dt + e f K{t)dt^ Pb{g) + e f K{t)dt 

Jo Jo Jo 

and ^ = 

Equation (2.2)3 can be reformulated in the following way being the modification of 
the entropy equation: 

- div f (1 + e^"') ^^^^ ^ Vs) + Kig)gv ■ Vs - Kig)v ■ Vg + div (v [' K(t)dt) 
(2 4:) ^ e^ / ^ Jo ^ 

^- (r.^ ^ \ S{v):Vv (l + e^"^)(e + e^),^ ,2 . ^ 
+ diY {KiQ)QV] ^ V- ^Vsr ml], 

with the "entropy" s defined as follows 

(2.5) s = ln^. 

The distinguished entropy will allow to control the positiveness of the temperature, what 
does not seem to be elementary working directly with equation of type (2.2)3. 
This system is completed by the boundary conditions at dQ, 

(1 + ^™)(e + ^)^ + - eo) + es = 0, 

(2.6) v-n = Q, Tk-{T{p,v)n) + fv-Tk = 0, k = l,2, 

1^ = °- 

on 

A key element in the limit passage from the approximative problem to the original 
one is the energy estimate giving information independent of the choice of function K, 
i.e. of the choice of the positive constant k — see (2.1): 

Lemma 1 Suppose solutions to (2.1)-(2.6) to he sufficiently smooth, i.e. g, v and 9 G 
W^iyt) for any q < 00, 9 > in fl. Let assumptions of Theorem 1 he satisfied. Then 

< g < k, / gdx < M and 
Jn 

(2-7) Mhhci) + \\Kig)g\\L,,in) + \\PiQ,9)\\L,in) + \\9\\L,m{n) + IIV^Hm^) 

+ / {e' + e-')da+\\Vs\\L,(n) <Ci\\F\\L^^n),M), 
Jon 

where the r.h.s. of (2.7) is independent of e and k, s = In^ and r — min{2, -^}. 
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Proof. The positiveness of the density and boundedness by k follow elementary from 
features of function K and the form of (2.2) i. The integration of this equation leads to 
the bound on the total mass. For details we refer to [11]. Let us prove the second part 
of (2.7) which is definitely more complicated. Multiply the approximative momentum 
equation (2.2)2 by v and integrate it over fl: 



(2.8) 



/ {2ij^B^{v) +Xdw'^v)dx+ / f\vQT\'^da+ / v-VPb{g)dx 
Jn Jon Jn 

= J K{g)gv ■ Fdx + J [ K{t)dt^ divvdx. 



To find a good form of the last term of the l.h.s. of (2.8) we use the approximative 
continuity equation (2.2)i. 



In' 7-1 

7 



v-VPb{g)dx^—^ / K{g)gv -Vg^-^dx 



- / [eA^ + ehK[g) — eg] g'^ ^dx 
1 Jn 

'-^ I [g-hK{g)]g^-^dx + e-^ I g^-^\V g\'^dx. 
~ ^ Jn Jn 

Thus the momentum equation gives the following inequality 

/ S{v):Vvdx+ I /Iv 0x1^^(7 + 67 / g'^''^\Vg\^dx + -^ I g^dx 
Jn Jan Jn 1 ~ ^ Jn 
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(2.9) 



-j (^j\{t)di^ediYvdx<c(l + j \K{g)gv F\dxy 



Integrating the energy equation (2.2)3 and employing the boundary condition (2.6) 1 we 
get 

(2.10) / {L{e)(e -eo) + es)da = [ (S(v) -.Vv- { [' K{t)dt)edivv)dx, 
Jan Jn^ Jo ^ 

since the integration by parts gives the following identity 

y K{g)gv ■ Ve - eK(g)v ■ Vg + div (v ^ K{t)dt^e 
+ diY{K{g)gv)e dx^ j j'^ K{t)dt^e diw vdx. 

Summing up (2.9) and (2.10) we get 

/ {L{9)e + es+)d(T + e-f [ g'^'^lX^ g\^dx + [ g'^dx 

(2.11) -^9^ . J^ . 7-Un 

< / es-da + C[l+ / \K{g)gv ■ F\dx\ 
Jan ^ Jn ' 

where s"*" and are the positive and negative parts of the entropy, respectively (s = 
S+-S-)- 
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We shall concentrate the attention on the first term of the r.h.s. of (2.11). Note that 
the control of the negative part of entropy s is not immediate. We integrate the entropy 
equation (2.4) over Q getting 

'm{9-9o) 

-I- fcii e uu -t- I i U] u IWUVU ■ ML 



(2.12) 



So 



(2.13) 



/ 



e 



/ 



e 



S{v):Vv (l + e"*)(e + 



e 



+ 
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\Vs\ 



dx. 



S{v) : Vv , (1 + ^"')(r + /?),^_^,2 



Vsr dx + 



[ K{g)gvV{s-]ng)dx< [ L{9)da + I 
Jq Jdn Ja 



es"^e ^ da. 



Let us look closer at the last term in the l.h.s. of (2.13). We have 
(2.14) 

— / K{g)gv ■ V(s — In g)dx = / K{g)gv ■ V In gdx — / K{g)gv ■ Vsdx = h + h, 
J^i Jo, Jci 

and employing (2.2) i we get 



(2.15) 



K{g)gv ■ V In gdx 



div{K{g)gv) In gdx 
\Vg\' 



/ [ — eAg + eg — ehK(g)^ In gdx = / (e ehK (g) In g + egln gjdx. 



The first term has a good sign, the second term has a good sign iox g < 1, too, and for 
^ > 1 is easily bounded by ehg. Similarly, the last term can be controlled by the term 
e g'dx. The proof was rather formal, as we do not know whether ^ > in Q. However, 
we may write K{g)v ■ V(^ + S) in (2.12) with 5 > and find an analogue of (2.15) with 
\n{g + S). Finally we pass with S ^ and get precisely the same information as above. 
Next 



(2.16) 



h — — K[g)gv-Vsdx— / [el^g — eg -\- ehK[g))sdx 
Jn J 

^ / { - eVgVs - eg\n9 + ehK{g)\n9)dx. 
Ja, 

Considering the r.h.s. of (2.16), we have 



(2.17) 



WgWsdx 



< e||V^||L2(n)||Vs||L2(Q) 



<lef / \^dx+ [ \Vg\'g^-'dx)+hvi 
4^ ^Jn Q Jn ' ^ 



\L2isiy 



Moreover, —egln9dx has a good sign for ^ < 1 and for ^ > 1 
(2.18) 

j^-eg{\n9ydx < e||^||L2(Q)||s+||L2(n) < |(||s+||Li(aQ) + ||Vs||l2(q)) + ^\\g'^\W{n) + C. 
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The last term of (2.16) can be treated as follows (one part has again a good sign) 
(2.19) / ehK{Q)\{\ne)-\dx <Ce [ \s~\dx<C + l I e\s~\e\'~\da s\\L^(n). 
Then combining (2.13) with inequality (2.11) and with (2.15)-(2.19) we obtain 

where 

H = C(^1 + J \K{g)gvF\dxy 
Thus from the growth conditions we deduce the following "homogeneous" estimates: 

To obtain a good information about integrability of the temperature we use the following 
Poincare type inequality 

( / l^^/T^^a; <C(0) ( / iV^^/^da; + ( / O'+^da 

\Jq ) ) \JdQ 

which can be proved elementary. Then the imbedding theorem leads to the bound 
(2.21) lje^"'dx] < i/V™ + ^i/a+i) 

To simphfy further calculations, we set Z + 1 = m. Note that we may allow also different 
values of /, however, for the prize that the further calculations become more technical 
which we try to avoid. 

We return to (2.9). Holder's inequality yields^ 



|2 

I I 'f I 

(2.22) 



^Il^i(n) + e7 / Q''-'\Vg\^dx+^ f Q^dx 



<C 1 + 



Jn J - ^ Jn 

[ \K{g)gv F\dx+ [ \9 [ K{t)dt\'^dx] . 
Jn Jq Jo / 



The next step of our estimation is the bound on Pb{g) which is necessary to estimate 
the r.h.s. of (2.22). We just repeat the method for the barotropic case, but here we shall 
obtain an extra term related to the temperature. 

Introduce # : Q — > defined as a solution to the following problem 

(2.23) t^v ''''' - Z £: {^'t-^)) - M /„ ^'C^)-^- 

^Note that we used Korn's inequality; for / = we therefore require that Q is not rotationally 
symmetric, for more details see [16]. 
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The basic theory to the stationary Stokes system gives the existence of a vector field 
satisfying (2.23) with the following estimate for a solution to (2.23) (for another possible 
proof, using directly estimates of special solutions to system (2.23), see [16]) 

(2.24) ||#|Ui(f,)<q|A||L.(f.). 

From the structure of Pb{g) and information that g^dx < M we easily get applying 
the interpolation inequality 

{Pb{g)} < S\\Pb{g)\\L,in)+C{5,M) for any 5 > 0. 

Multiplying the momentum equation (2.2)2 by employing (2.22) and (2.24), we con- 
clude after standard estimates of the r.h.s to (2.2)2 

(2.25) I in(^) I <c{l + \K{g)gv ® v\^dx + jj^ K{t)dt\^dx^ . 
As 

(2.26) ||n(^)|lL(n) > C i^jj^K{g)gf^dx-r (^j'^ K{t)dt^ , 

recalling that 27 > 6, we get a bound for the first integral in the r.h.s. of (2.25) 
(2.27) 

/ \K{g)gv ® v\^dx < c\ \v\\]J,^^^ \\K{g)g\ |i,(n) 
Jo. 

2(7-3) 107 6(27-1) 

< cii^ii^,i(o)iii^(^)^iij7nTii^(^>)^iiSo) ^ ^\\p^mim+c^^^M)\Mmi:i ■ 

Hence a suitable choice of 5 in (2.27) simplifies (2.25) to 

/ 6(27-1) r \ 

(2.28) ||n(^?)|lL(o) <C[1 + ||t;||^^i74 + \e K{t)dt\'dxj . 

The last integral can be viewed by (2.26) in the form 
(2.29) 

/ 3 2^-1 , r rQ ^ 



/ 

Jo 



K{t)dt\\L,,in) + ||i^(^)^||L,,(n) <c(^l + \\v\\lt^^^; + JjO j^' K{t)dt\'dx 



Within our estimation we concentrate on a precise specification of powers of norms. Then, 
due to our growth conditions we shall be able to construct the desired bound (2.7). 
The last integral in (2.29) can be treated as follows (we need rn > ^ and m > -^^zrj) 

Q 



(2.30) 



1 1/7 

3m-2 



ii^£x(t)dtii^5^)<ii^ii^(:(^)ii£x(t)dtii 

1 re (3m-2)7-3m pQ 3m + 2 

^ll^llL(n)ll I mdt\\,Ti^^^-'' II I ^W^^IIlSni" 

so (2.29) and (2.30) with the Holder inequality imply 

re , 3 27-1 3m 27-1 ^ 

II i^(t)rft||wn) + ||i^(^)^|lwn) <C'(l + ||^||J?-; + |^ 
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(2.31) II / K{t)dt\\r.,^^n) + ||i^(^)^||L.,(n) < C (l + \\v\\llj^^ + H^''''^ 



Applying the inequality for the temperature — (2.21) — we obtain (recall that we put 
l + l = m) 

3 27-1 3 2-,-l 

37-4 

We have to estimate H; it holds 

\K{g)gvF\dx < \ \v\\Lein)\\K{g)g\\L,/,in)\\F\\L^^n)- 
Using the interpolation between 1 and 27 as above leads to the following bound 

(2.32) [ \K{g)gvF\dx<C{M)\\v\\HHn)\\K{g)g\\^]. 

Inserting this inequality to the r.h.s. of (2.31), recalling that m > ^ and applying the 
standard Holder inequality we obtain from (2.31) estimate on the density 

(2.33) \\ J' K{t)dt\\,,^^n) + ||i^(^)^IU.,(n) < c(l + ||t;|||f|^ + H^lli^f^). 

As we can see later, the first term is the most restrictive. So by (2.32) and (2.33) we 
conclude (for m > |^E|) 

(2.34) f \K{g)gvF\dx <c(l + \\v\\^^^y 
Hence we obtain from (2.21) 

(2.35) ||^||L3.(n)<c(^l + ||7;|||,|f 
From (2.30) we easily see that 

3m + 2 7 



(2.36) 11^ \ mdt\\^,^^^<C\n^,^^^^\\ I K{t)dt\\,l-^^''-' 



Summing up inequalities (2.22), (2.34) and (2.36) we obtain the main bound on the 
norm of the velocity 

37-3 2 37-3 , 2 3m+2 

IdIP, < r" I 1 -I- I IdI P^"* _|_ I U,! I 37-4"^m 37-4 



|f/i(a) ^J- -r Mt^Miji(f2) 11*^11^1(0) 

The above bound implies the a priori bound 
(2-37) ||t;||i/i(Q)<C(||F|U^,M), 

provided suitable dependence between 7 and m holds, which can be described by the 
sufficient condition (7 > 3) 

(2.38) m > 

Note that as wc take 7 near 3 then m > 4 and for 7 = 4 we have m > Moreover, 
the above needed conditions m > |^Ei) ^ > | m > 3(^zij clearly less restrictive 
than (2.38). 

Bound (2.37) implies immediately the a priori estimate (2.7), since it follows from 
(2.20) with (1.11), (2.28), (2.33)-(2.36), together with (3.7). □ 
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3 Existence for the approximative system 

The aim of this section is to show that for any e > and k > there is a solution to 
the approximative system (2.2)-(2.6). We prove 

Theorem 2 Let the assvmptions of Theorem 1 he satisfied. Moreover, let e > and 
k > 0. Then there exists a strong solution {q,v,s) to (2.2) such that 

Q e (O), V e Wl{Q) and s e for 1 < p < oo. 

Moreover < g < k in fl, gdx < M and 

(3-1) \MwU^) + V~e\\Vg\\L,in) + ||V^|U,(n) + ||^||L3^(n) < C{k), 



e. 



where 6 = e^ , r = min{2, -^^} and the r.h.s. of (3.1) is independent of the parameter 

The proof of the existence to the approximative system (2.2) will follow from the 
standard application of the Leray-Schauder fixed point theorem. It will be split into 
several lemmas. First we consider the continuity equation. We denote for p e [1, oo] 

Mp^{we W^{ny, w-n = at dn}. 

We have 

Lemma 2 Let q > 3. Then the operator 

S:Mq^ W^(n) for 1 < p < oo 
such that S{v) = g, where g is the solution to the following problem 

eg — e/S.g — ehK{g) — div(ir(^)^v) in Q, 

(3-2) dg 

Tr^ = at dn 

on 

is a well defined continuous compact operator from Mg to Wp{fl), 1 < p < oo. In 
particular, the solution to (3.2) is unique. Moreover 

(3.3) MWim < C{k, e)(||t;||^(-i(^) + 1), / = 1, 2. 

Proof. It follows from [11], Proposition 3.1 (there, the two dimensional case was 
considered). See also [16]. □ 
Next, we define the operator 

T:MpX W^(n) ^MpX W^(n) such that T(v, s) = (w, z), 
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where {w, z) is the solution to the following system 
(3.4) 

- div S(ti;) = -- diy{K{g)gv <^ v) - -K{q)qv ■ Vv - VP{g, e*) + K{q)qF 

-div((l + e'"')(e + e')V;2) = ^{v) : Vt; - div [v [\{t)dt)e' 

^ Jo ^ 
- div {K{g)QV)e'' - e'K{Q)QV ■ Vs + e'K{Q)v ■ Vg ^ 

w ■ n — 0, n ■ S{w) ■ Ti + fw • T; = for Z = 1, 2 

(1 + e"^^)(e + e')Vz + = -L(e")(e" - ^o) 



> in 



at dVt, 



where q = S{w) is given by Lemma 2. 

Our aim is to apply the Leray-Schauder fixed point theorem. Thus we need to verify 
that T is a continuous and compact mapping from Mp x Wp{Q.) to Mp x Wp{^) and that 
all solutions satisfying 

(3.5) tT{w, z) = {w,z), te [0, 1] are bounded in Mp x W^{n). 
First we easily have 

Lemma 3 Let p > 3 and all assumptions of Theorem 2 be satisfied. Then T is a contin- 
uous and compact operator from Mp x Wp(Q) to Mp x Wp{D,). 

Proof Note that for e > the system (3.4) is strictly elliptic. Since p > 3, the 
VTp (f2) -space is algebra, thus the r.h.s. of (3.4) belongs to the Lp-space (the boundary 

term belongs to Wp'~^^'^{dQ)). The coefficients in the operator in the l.h.s. of (3.4)2 are of 
the C^"'""(n)-class. Hence the standard theory for elliptic systems gives us the existence 
of the solution to (3.4) in Mp x Wp{il) with the following bound 

||w^||w|(a) + < C'(l|e1lci+«(n)) (lithe r.h.s. of (3.4)i| 

+||the r.h.s. of (3.4)2| |Lp(n) + ||the r.h.s. of (3.4)4| l^^i-i/p^g^^)) 

which guarantees us the uniqueness and the continuous dependence on the data. Moreover 
the r.h.s. of (3.4) is at most of the first order of sought functions. Thus this structure 
implies the compactness for the map T. □ 
Next we consider a priori bounds for solutions to (3.5). 

Lemma 4 All solutions to problem (3.5) in the class Mp x Wp{p) satisfy the following 
bounds 

(3.6) Q<Q<k, |k||//i(n) + ||^||L3^(n) + ||V^||L,(n) + V~e\\WQ\\L,(a) < C{k), 

where r — min{^|^, 2}, 6 — and the constant C(k) is independent of e and t e [0, 1]. 

Proof. We may basically repeat estimates of Lemma 1 from the previous section. 
However, on the one hand, we are in a simpler situation as we can use bounds which 
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depend on k, i.e. on the bound of the density (they may be proved analogously as in 
[11]), on the other hand, we must control the behavior of all norms with respect to t. 

Thus, repeating steps (2.8)-(2.13) for the case t — 1 (the corresponding terms are only 
multiplied by t) we finally get 

(l-t) f S{w):Vwdx+ f f{wQTfda+ f + ^")(^ + ^) |y^|2^^ 
Jq Jan Jn ^ 

+e / [z+(l-e-^+) + |^_|(el^-l-l)]da + i / [L{e)e - L{e)e^ + -^^^ - L{e)\da 
Jan Jen ^ 

<t J (^K{g)gw-Vz-K{g)w-Vg^dx + tc(l + J \K{g)gw ■ F\dx^ , 
where g = S{w). 

We may now repeat the arguments between (2.14)-(2.20) (all the corresponding terms 
are only multiplied by t) and we finally get 



Jn f Jn Jan^ ^ 

\K(g)gvF\dxy 



<tC 1 



As < g < k,wc easily get (the Poincarc inequality is just the same as in the previous 
section), after dividing by t (the case t = is clear; recall also m = I + 1) 

||^||L3.(n)<C(l + ||w;|U,(n))^/"^ 
and from an analogue to (2.22) also 

lkllli(Q)<c(i + 11^11 

As m > 1, it implies 

IHIhhq) + \\0\\L,^{n) < C{k). 

Further, if m > 2 then due to the control of and \W0\0^^ in ^2(^2) we have also 
bounded in the same space. For 1 < m < 2, 

(3.7) IIV^IU 3. (Q) < ll|V^|^^|U.(n)||^ILl(n)- 

m+1 

Finally, multiplying the approximative continuity equation by g and integrating by parts 
we get 

e / (|V^|^ + ^^)cix < e / hK{g)gdx+ 1(1 K {t)tdt)\ diY w\dx, 
Jn Jn Jn^Jo ' 

from where we deduce the bound for ^/e|| Vf)||L2(n)- n 
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To conclude, we verify the bound on {w, z) in W^iVt) x W^{Vt), p < oo, independently 
of t. We apply the bootstrap method to system 

— div S{w) —t — div{K{g)gw (S>w) — -K{g)gw ■ Vw 



(3.8) 



-VP{g,e^) + K{g)gF 

■dw{{l + e"''){e + e')Vz) S(mj) : Vw; - div (w j' K{t)dty 
- div {K{g)gw)e' - e^K{g)gw ■ Wz + e^K{g)w ■ V^] 



> in Q,, 



w ■ n — 0, n- S{w) ■ ti + fw ■ ti — Q for Z = 1, 2 

(1 + e"'^)(e + e')Vz + ez = -tL{e'){e' - 60) 



at dil, 



where g = S{w) given by Lemma 2. Note first that due to bounds from Lemma 4 we 
have 

ll^llvyicn) < C 

as K{g)gw ® w is bounded in 1/3(0). Thus w is bounded in any I/g(f2), q < 00 and 
the most restrictive term is VP(^, e^). As e^ = ^ is bounded in L3^(Q), g in Loo(il), we 
deduce the bound 

< C and consequently also ^ C*. 

Note that the constant in the estimate for w is independent of e. 
Next, we rewrite equation (8.8)2 as follows 



(3.9) 

with 
(3.10) 



-A^z) = t 
— div I w 



S{w) : Vw + e^K{g)gw ■ Vz - e^K{g)w ■ Vg 



K{t)dt^e' - diY {K{g)gwy 



in Q, 



dn 



^ -ez - tL{e^){e^ - Oo) at dQ 



^z) = / (l + e"*^)(e + e^)dT. 
Jo 



We multiply (3.9) 1 by $ and integrate over fl. It leads to 



IV$| 



L2{n) 



[ {tL{e'){e' -9o)^ + ez^)da<C\\theT.h.s. of (3.9)i| |i |$| 
Jan 



Le{n)- 



It is not difficult to realize that the most restrictive term on the r.h.s is e^K{g)gw ■ Vz e 

L 3m (fi), where > | for m > 1. 

Let us look at the boundary terms. Note that ^{s) ~ es for s — > —00 and $ ~ g(m+i)« 
for s — > +00. Thus 

/ [tL(e^)(e^ - eo)^ + es$]7{$<o}da > P||L.(an) - C 
Jan 
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and 

/ [tL{e'){e' - Oo)^ + es^]l{^>o}da > \MMdn) - C. 
Jaa, 

Thus, the estimates above yield ^ ^ with C independent of t which implies 

IK^'heici) = ||e("^+^)^||Ls(f2) < C and also ||V^||L,(n) = ||e^V;^||L,(o) < C. 

Now, it is not difficult to verify that from (3.9) we get ||^||vK2»(n) < C with p* — min{^, 2} 
(thus e^Vz e L2{0,) and Vty e L3^(0)). In particular, 

ML^m + WOh^in) < C, WVzh.in) + ||V^^|U,(o) < c 

for 1 < q < q* = jzp^ > 3. Thus from the approximative momentum equation we get 
(V(^6') e Lq*{fl)) the bound ||iy||iy2^(f2) ^ C and from the energy/entropy equation also 

The imbedding theorem yields ||V^||L^(n) + ||V^||L^(n) < C which finally gives as above 

||w'||w2(n) + Ikl|w2(n) + ||6'||w2(n) < C, 1 < r < oo 
with C independent of t. This finishes the proof of Theorem 2. 

4 Effective viscous flux 

In this part we investigate the properties of the effective viscous flux. Estimates (3.1) 
from Theorem 2 guarantee us existence of a subsequence e — > 0"^ such that 

Ve^v in Loo(fi), 
Qe ^* Q in Loo(fi), 

Ph{Qe) ^* Pbjg) in L^{n), 
, , ^(gc)gc ^* K{g) g in L^{n), 

^ ■ ' K{g,) K{g) in Lo,{n), 



j^^ K(t)dt K{t)dt in L^{n), 

e,^e in W^{n) with r = min{2, ^}, 
9e ^ in Lq{0,) for q < 3m. 



Passing to the limit in the weak formulation of our problem we get 



(4.2) div{K{g)gv) = 0, 



(4.3) K{g)gv ■ Vv - div {2nD{v) + i/(divy)/ - Pb{g)I - e{j'^ K{t)di)I^ = K{g)gF, 
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(4.4) - div((l + e"')Ve) + e{dw v r K{t)dt) + div(K{Q)Qev) = 2/i|D(v) ^ + i/(div v)^ 

Jo 

together with the boundary conditions (1.9)-(1.10). Recall that (4.2)-(4.4) is satisfied in 
the weak sense, similar to Definition 1. 

In what follows we must carefully study the dependence of the a priori bounds on k. 
We have 

Lemma 5 Under the assumptions of Theorems 1 and 2, we have 

(4-5) WQeWhooiO) < k and \\ve\\wl^{Q.) < C{1 + k^^^). 

Proof. The bound on the density follows directly from Theorem 2. We therefore 
estimate the velocity. If we write (2.2)2 in the form 

-div5(i;) = -v(Pk{Qe)+e,{ r K{t)dt))+K{g,)g,F 

1 \ 
-- diY[K{ge)g^v^ v^] - -K{g^)g^v^ ■ Vv^, 

we immediately see that 

llv.llwi^Cn) < C{\\K{g,)g,v, ® ■y,||L3^(n) + \\K{g,)g,v, ■ Vv,\\l 3^ (q) 

m+1 

+ ||n(^e)||L3^(a) + ||^.( / K{t)dt)U,^^n) + \\K{g,)g,FU 

Jo "»TT 

Note that due to the bound of the temperature we cannot expect e-independent estimate 
for q > 3m. The bounds on the density and temperature yield 
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||n(^e)||L3.(n) < \mQe)\\l:^a)\\Pb{Qe)\\C(n) < Ck^^, 

while 

11^, ( r K{t)dt)U,^^n)<Ck. 
Jo 

Note that for m and 7 satisfying assumptions of Theorem 1, 7 ^™^^ > 1. It remains to 
estimate the convective terms (CT.) 

L 3m 

< C'll^6|Uoc(n)(||^^6||i6^(f2) + IIVV^IIl 3^ (Q)||^^6||Loo(n)) 

m+1 

for m > 2, while for m < 2 the last term is replaced by || VVell 1,2(0) H'i^ellL 6m {^)- Using 

2 — 771 

the fact that for 6 < g < 00 

||^.|U,(n) < C||t;.||2,(n)ll^e||^3^(^) with - = i + (1 - «)(^ - 3) 

and for 2 < r < 3m 

||V-y,|U,(n) < \\vTL,i^)\\^v,\\\-Zin) with ^ = f + 
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we end up with 



C.T. < C\\ge\\L^(n)\\Ve\\^T[^-j\\Ve\\^Tj^y 

Note that < 1. Thus we may use the bound on and Young's inequahty yields 

, 7 3m-2 , 3m-2 1 

< C(l+/C3 m )+Ck m ^-\\v^\\^^^^^y 

As 7 > 3, the lemma is proved. □ 
Before using the above proved bounds, we show one useful result which in particular 
implies that the limit temperature is positive. 

Lemma 6 There exists a subsequence {s^} such that 

^ s in iv2(il), 

subsequently, 

9 in Lq{Q), q < 3m with 9 > a.e. in Q. 
Proof. Recall that from the energy bound we have the following information 

[ \Vse\'^dx+ I {e'^ + e-'^)da < C 
Jn JdQ 

which in particular gives 

/ \Vse\^dx+ / S^d(7 < C. 
Jn Jdn 

Thus we are allowed to choose a subsequence Sg — > s in L2{fl). Recall also that 9^ = e^^ 
and 9^ ^ 9 strongly in Lr{fl), r < 3m. Hence by Vitali's theorem (for a subsequence, if 
necessary) 

e^' in Lr{fl) and 9 = e'' with s e L2{fl). 

Thus ^ > a.e. in Q as s > — oo a.e. in Jl. □ 

A crucial role in the proof of the strong convergence of the density is played by a quan- 
tity called the effective viscous flux. To define it, we need the Helmholtz decomposition 
of the velocity 

(4.6) V = V0 + rot A, 

where the divergence-free part of the velocity is given as a solution to the following elliptic 
problem 

rot rot A = rot v = uj in fl, 

(4.7) div rot A = in fl, 

rot A • n = at dQ. 



18 



The potential part of the velocity is given by the solution to 

A6 = divv in /" , , 

(^•8) 1 = atafi, //''^ = °- 

The classical theory for elliptic equations gives us for 1 < g < oo 

llVrot < C||w||L^(n), HV^rot A||i^(n)|| < C\\uj\\w^^n), 

||V2(/)|U,(n) < C||divv||L^(n), ||V=^(/)|U,(n) < C|| divy||w^i(n). 

The properties of the slip boundary condition enables us to state the following problem 

-^AUJ, = rot {K{Qe)QeF - K{Qe)QeVe ' Vv^ 

'-]^thK{Q^)Ve + ]^tQeVe) - mti^e/^Q^v^) := Hi + H2 in Q, 

a;,-Ti = -(2x2-///^K-T2 ataO, 
We • T2 = (2xi - • Ti at aQ, 

divWe = at 9Q, 

where Xfc are curvatures related with directions Tfc. For the proof of relations (4.9)2,3 ~ 
see [10] or [12]. 

The structure of oJg gives us a hint to consider it as a sum of three components 
(4.10) u}, = u}l + uj\ + u}l, 

where they are determined by the following systems 

-/xAuJ° = 0, -^Aa;i = H 

a 11^ ■ Tl = -(2X2 - //^)«. ■ T2, u\-Ti = Q, 

^^■''> a;0 ■ T2 = (2x1 - ///x)t^e • Ti, u:\ • = 0, 
diva;° = 0, divw.^^O, 

Lemma 7 For the vorticity uj^ written in the form (4-10) we haveP' 

||u;?|lMn) <C(A;)eV2 /or 1 < r < 2, 

(4.12) 4 2 

+ \\<^\\\w^,{U) < C(l + A;^+^(3-5)) for2<q< 3m. 

Proof. First, let us consider 0;°. Take ao any divergence-free extension of the 
boundary data to Wg, e.g. in the form of a solution to the following Stokes problem 

—/iAcxq + Vpo = in 
div oto — in Jl, 

(4.13) Cto-Ti = -{2X2- f/l^)Ve-T2 Bit dQ, 

ao • T2 = (2xi - f/l^)v, ■ T2 at dQ, 
cxq - n— at dQ. 



. n 


= J/2 


in r2, 


= 


at d^l, 


■ T2 


= 


at (9^2, 


div u?g 


= 


at dQ. 



^Note that wc can prove that = o(e) for e — ^ 0+ for any r < 3m. As wc do not need it and 

the proof of the rate is shghtly more comphcated, we skip it. Analogously we may consider the other 
inequality also for g < 2, with different powers of k. 
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Note that e W^,n^^^^"'\dfl), thus (Xq e W^^{fl) with the estimate 

||tto|Ui{n) < C||ve||wi(n), l<q<3m. 
Thus we may transform the system for a;^ to the form 



(4-14) . ,0 



—fj,A{u}'^ — ao) = /xAcKo in Q, 
(u;° - ao) ■ Ti = at dQ, 
(a;° - ao) • T2 = at 91^, 
div(u;^' - ao) = at dQ. 



Note that Ackq e W3^(0). Here ^{fl) denotes the dual space to 

{/ G w^iQ) n {/ ■ n = / ■ T2 = at an}}. 

As the system for c«;° has the same structure as that for u)], we get (see [18], [19]): 
^ C\\^^i\\w,-\n) and ||u>°||w'i(f2) < C\\v,\\wi{n), l<q<Sm. 

Analyzing the form of Hi we see that the only not elementary term is the convective 
one; so we obtain 

I \<^l\\w^{n) < C(l + \\K{g,)g,v, ■ Vv,| |L^(n)). 
We easily see that for g > 2 

\\K{ge)geVe ' Vt;,||L^(n) < k\\Ve\\L^{n)\\^Ve\\Lg{n)- 

Using interpolation inequalities as in Lemma 5 we prove that 

2(m— 1) m 6m — 2q 3m(q — 2) 

\\K{g,)g,v,-Vv,\\,,ia) < ^^^ll^.llj^ II Vt;,||S^) || Vt;,||g^ || V^;.||^ 

Evidently, the estimate for u}^ is less restrictive. 
Similarly, for uj^ we have 

\\uj^,\\L,{n) <C\\eAg,v,\\^-i^^^<Cesnp\ / Ag,v,(j)dx\, 

where the sup is taken over all functions belonging to Wq{Q) with 1/p + 1/q' = 1. 
Prom the continuity equation we know that 

\/^l|Vf?,|U,(f,) < c{k). 

(For q > 2 we have only e|| V^g||z,g(n) < C.) As g < 2, 

||w^||l,(q) < C'e(||V^e||L2(n) 11^,111,^(0) + ||V^>e||L2(n)||Vv,|| ) < C{k)e-^. 
The lemma is proved. □ 
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We now introduce the fundamental quantity — the effective viscous flux — which is 
in fact the potential part of the momentum equation. Using the Helmholtz decomposition 
in the approximative momentum equation we have 

V(-(2/i + u)A(f), + P{g„ e,)) = /xArot A, + K{g,)g,F 

1 11 
-K{ge)geVe • Vve - -ehK{ge)ve + -egeVe - -eAg^Ve- 

We define 

(4.15) = -{2ii + u)A(j), + P{g,, 9,) = -(2/i + u) div + P{g,, 9,) 
and its limit version 

(4.16) G = -(2/i + z/) diYV + P{g,9). 

Note that wc arc able to control integrals j^G^dx = j^P{ge,9e)dx and J^Gdx = 
P{g, 9)dx, where P{g,9) = Pb{g) + 9 {J^ K{t)dt) . 

The result of the lemma below gives the most important properties of the effective 
viscous flux, guaranteeing the compactness of {Ge} as well as the pointwise bound of the 
limit in term of the parameter k from definition (2.1). 

Lemma 8 We have, up to a subsequence e ^ O"^.' 

(4.17) G^^G strongly m L2{^) 
and 

(4.18) \\G\\L^<C{7]){l + k^+'^"'+'^) for any 7] >0. 
Proof. The function G^ can be naturally decomposed as 

G,^Gl + Gl 

where G^dx — and VG^ = —\eAg^v^ — //rot u?^. Thus 

\\GI\\l,{Q) < C(e||A^,'y,||^^-i(j^) +/i||rot w^H^-i^^^). 
Using Lemma 7 we see that 

\\Gl\\L,in)<C{k)el, l<q<2. 

Next, using again Lemma 7 and calculations in its proof, we immediately see that (recall 
that I J^G4x\ < C) 

(4.19) ||G',i||w^i(n) <C(l + A;^+^(^-|)) for 2 < g < 3m. 
Thus we have, at least for a subsequence 

G\ G^ in L~(0) and G^ ^ in L^i^). 

Therefore 

G, = G\ + ^ in L^iSl), 1 < g < 2 

and due to the definition, G^ = G. Finally, choosing g = 3 + f/ in (4.19) 

||G|U^(n) < G(g)||G|Ui(n) < G(g) sup ||G,iH^i(n) < G{r)){l + k'^h^"^) 

e>0 

with 77 > 0, arbitrarily small if 77 is so. This finishes the proof of Lemma 8. □ 
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5 Limit passage 

In this section we apply the properties of the effective viscous flux shown in the previous 
part. First we prove a result characterizing the sequence of approximative densities. 

Theorem 3 There exits a sufficiently large number ko > such that for k > ko 

(5.1) ^(^-3r-||G|U^(n)>l 

and for a subsequence e ^ 0^ it holds 

(5.2) lim \{xen: g,{x) >k-3}\^0. 



In particular it follows: K{q)q = g a.e. in ^l. 

Proof. We define a smooth function M : Rq — > [0, 1] such that 

( 1 for t<k-3 

M{t) = <^ e [0, 1] for k-?,<t<k-2 
[ for k-2<t 

and M'{t) < for t e {k - 3,k - 2). 

We follow the method introduced in [11]. First we multiply the approximative conti- 
nuity equation (2.2) i by M\q^) for I e N getting 

r ( rer\x) \ 
j \j tlM^-\t)M'{t)dt\ diyv, > R, 

with — > as e — > 0, as 

e [ M^{g,)Ag,dx^ -el [ M^-\g,)M'{g,)\V g^l^dx > 0. 
Next, recalling definitions of and M, we obtain 

-(fc-3) J (^£''^\M^-\t)M'{t)dt^P{g,,e,)dx 



< k 




lM^~\t)M'{t)dt^G,dx 



+ Re. 



'n ^ Jo 

Thus the properties of M lead us to the following inequality 

^ f {l-M\g,))P{g,,e,)dx< f {l-M\g,))\G,\dx+\R, 

1^ J{ge>k-3} J{ee>k-3} 

From the explicit form of the pressure function (2.3) we find 

^(^-3)^|{^. > A;-3}|-^||P(^„^e)||L.(Q)||M'(^,)||i,(f,) 
< \\G\\L^in)\{ge > A; - 3}| + \\G - ^11^,(0) + \Re\- 
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But by Lemma 8 - the inequality (4.18) - we are able to choose ko so large that for 
all k> kowe have (5.1), since 7 > 3 and ||G'||l^(o) < Cr,{l + k^+h+^) with 0<ri<^. 
Hence we get 

(5.3) \{xen: g,{x) > k - 3}\ < C {\\M\g,)\\L,i{,,>k-3}) + \\G - G',||i,(n) + \Re\) ■ 
Now, let us fix 5 > 0. Then there exists eo > such that for e < eo 

(5.4) C{\\G-G,\\L,ia) + \Re\)<S/2. 

Having e fixed, we consider the sequence {M'(^g)7{g^>fe_3}};eN, where Ia is the character- 
istic function of a set A. We see that it monotoncly pointwise converges to zero. Thus by 
the Lebesgue theorem we are able to find I = l{e, 5) such that 

(5.5) C\\M^(6e)\\L,i{,.>k-3})<5/2. 
From (5.3), (5.4) and (5.5) we obtain 

(5.6) \im\{x e Q; Qeix) > k - 3}\ < 5. 

As 5 > can be chosen arbitrarily small, Theorem 3 is proved. □ 
Thanks to Theorem 3 we are prepared to present the main part of the proof, i.e. the 
pointwise convergence of the density. 

Lemma 9 We have 

(5.7) / P{g,9)Qdx < / Ggdx and / P{g,9)Qdx= / Ggdx 
Jq Ja Jn Jci 



consequently, P{g, 0)g = P{g, 9)g and up to a subsequence e — > 0"^ 

(5.8) Qe ^ Q strongly in Lq{fl) for any q < 00. 

Proof. Due to Theorem 3 we are able to omit K(g) in the limit equation. For details 
we refer to [11] - section 4, consideration for (4.16). 

Examine the approximative continuity equation (2.2) 1. We use as test function ln(^g-|- 
S) and passing with S ^ we obtain 

(5.9) / K{g^)v, ■ Vgedx > eC{k), 

Jq 

thus Theorem 3 implies 

(5.10) — / gedivVedx > Re- 

Jn 

Applying (4.15) to (5.10), passing with e ^ 0, then by the strong convergence of G^ - see 
(4.17) - we conclude that Gg = Gg, so the first relation in (5.7) is proved. 
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Next we consider the limit to the continuity equation, i.e. div(^'u) = 0. Testing it by 
In Q with an apphcation of Priedrich's lemma to have possibility to use test functions with 
lower regularity we obtain (for details see [11]) 

QdlYvdx — 0. 

The definition of G - (4.16) - shows the second part of (5.7). 

Due to elementary properties of weak limits we get gP{g, 9) < P{g, 6)g a.e. in Q, but 
(5.7) implies J^{P{g, 9)g — P{g, 9) g)dx < 0, hence 

gP(g,9) = P(g,9)g a.e., i.e. g'y+^ +^0 ^^g + g^e a.e. 

However, ^ g'^g and g'^9 > g'^9, so 

^7+1 = '^g a.e. and g'^9 = g^9 a.e. 

By Lemma 6 the temperature ^ > a.e., we conclude g"^ — g^ and for a suitably taken 
subsequence 

(5.11) lini||^,-^||i^=:^-^2^0. 

Thus the limit (5.11) implies ge —>■ g strongly in L2(fi) and by the pointwise boundedness 
of ge and g we conclude (5.8). □ 

Next, we would like to study the limit of the energy equation. The first observation 
concerns the velocity, we obtain the strong convergence of its gradient. 

Recall that from Theorem 3 and due to the strong convergence of the temperature it 
follows 

P{ge, 9e) p(g, 9) strongly in L2(ll), 

hence (4.17) implies 

(5.12) diWe^divv strongly in L2(0). 
Additionally we already proved that 

(5.13) rot t;^ — >^ rot v strongly in ^2(^2), 

since we observed that the vorticity can be written as sum of two parts, one bounded in 
Wq{VL) and the other one going strongly to zero in L2(J1). 

The regularity of systems (4.7) and (4.8) and convergences (5.12) and (5.13) imply 
immediately that 

— > V strongly in H^{Vt). 

In particular, we get 

(5.14) S{ve) : Vve S{v) : Vv strongly at least in Li(Q). 



If, 
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This fact will be crucial in considerations for the limit of the energy equation. Recall 
that 



(5.15) 



^ g in Lg{fl) for q < oo, 
— > V in Wg{Q) for q < 3m, 
9f 9 in Lg{Q) for q < 3m, 
^e-OinW^^^^^m. 



Consider the weak form of (2.2)3. For a smooth function we have 

e + 9. 



(5.16) 



/ (1 + 9^f-^V9, ■ V(t>dx + / L{9,){9, - 9o)(l)da 
Jn Jan 

J K{t)dtjv, 



V{9,(l>) + K(g,)g,v,-V(9,(l>) 



dx 



+ 



/ K{g,)g,v, -^9,(1) + div {9, V, (I)) / K{t)dt dx ^ S{v,) : Vv,(j)dx. 
Jn Jo Jn 



Thanks to (5.15), 



[i + e: 



e + 9. 



V9e^ {1 + 9"')V9 inLi(f]). 



(5.17) 



Passing to the hmit with the last four terms of the l.h.s. of (5.16) we get 

[-gvV{9(p) - gvV{94>) + g4>vV9 + diy {94>v)g] dx 

— / [—g9v ■ W(f) + g9 div V(f)\ dx. 
Jn 

In (5.17) we essentially used the strong convergence of the density. 

To control the behavior of the boundary term we note that due to (5.15)4 we see that 
9\an strongly in Li^i{dfl). Thus recalling (5.14) we get at the limit 



(5.18) 



/" (1 + 9"')V9 ■ V(j)dx + [ L{9){9 - 9o)da - [ g9v ■ Vcpdx 
Jn Jan Jn 

— / S{v) : Vv(j)dx— I g9 divv(f)dx. 
Jn Jn 



To conclude, note that we may show that the limit functions 9 and v belong to Wp{fl) 

for any p < 00. To see this, we introduce the function ^(9) = Jq{1 + t"^)dt, similarly as 
in Section 3, formula (3.10). Thus from (5.18) we immediately see that 9 G Loo(^^) and 
V e Wp{fl) for any p < 00. Using this fact once more in the energy equation, we observe 
that 9 e Wp{il), p < 00. Theorem 1 is proved. 
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